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Abstract. We study a simple continuous-time multi-agent system related to Krause's model of 
opinion dynamics: each agent holds a real value, and this value is continuously attracted by every 
other value differing from it by less than 1, with an intensity proportional to the difference. 

We prove convergence to a set of clusters, with the agents in each cluster sharing a common 
value, and provide a lower bound on the distance between clusters at a stable equilibrium, under a 
suitable notion of multi-agent system stability. 

To better understand the behavior of the system for a large number of agents, we introduce 
a variant involving a continuum of agents. We prove, under some conditions, the existence of a 
solution to the system dynamics, convergence to clusters, and a non-trivial lower bound on the 
distance between clusters. Finally, we establish that the continuum model accurately represents the 
asymptotic behavior of a system with a finite but large number of agents. 

1. Introduction. We study a continuous-time multi-agent model: each of n 
agents, labeled 1, . . . , n, maintains a real number ("opinion") Xi(t), which is a contin- 
uous function of time and evolves according to the integral equation version of 

±i{t)= J2 (Xj(t) - Xi(t)) . (1.1) 

j: \xi(t)-x } (t)\<l 

This model has an interpretation in terms of opinion dynamics: an agent considers 
another agent to be a neighbor if their opinions differ by less than 1, and agent opin- 
ions are continuously attracted by their neighbors' opinions. Numerical simulations 
show that the system converges to clusters inside which all agents share a common 
value. Different clusters lie at a distance of at least 1 from each other, and often ap- 
proximately 2, as shown in Figure ITTTI We focus on understanding these convergence 
properties and the structure of the set of clusters, including the asymptotic behavior 
for large n. 

Observe that the agent interaction topology in (|2.ip explicitly depends on the 
agent states, as Xj(t) influences x.^t + 1) only if |xj(f) — Xj(t)\ < 1. Many multi- 
agent systems involve a changing interaction topology; see e.g. [1, 10, 11, 16, 19, 20], 
and [17, 18] for surveys. In some cases, the interaction topology evolves randomly or 
according to some exogenous scheme, but in other cases it is modeled as a function 
of the agent states. With some exceptions [6,7,12], however, this state-dependence is 
not taken into account in the analysis, probably due to the technical difficulties that 
it presents. 

To address this issue, we have recently analyzed [3] one of the simplest discrete- 
time multi-agent systems with state-dependent interaction topologies, namely, Krause's 
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Figure 1.1. Evolution with time of the values Xi(t) for 1000 agents, with initial values randomly 
and uniformly distributed on [0, 10] . Observe the convergence to 4 clusters separated by slightly more 
than 2. 

models of opinion dynamics [13]: n agents maintain real numbers ("opinions") Xi(t), 
i = 1, . . . , n, and synchronously update them as follows: 

^ {t+1) = ^5W!<l^ i 

2s j: \xi(t)- Xj (t)\<l 1 

This model was particularly appealing due to its simple formulation, and due to some 
peculiar behaviors that it exhibits, which cannot be explained without taking into 
account the explicit dynamics of the interaction topology. Indeed, a first analysis using 
results on infinite inhomogeneous matrix products, as in [10,14], shows convergence 
to clusters in which all agents share the same opinion, and that the distance between 
any two clusters is at least 1. Numerical simulations, however, show a qualitative 
behavior similar to the one shown in Figure [TTTI for the model (|2.1[) : the distance 
between consecutive clusters is usually significantly larger than 1, and typically close 
to 2 when the number of agents is sufficiently large, a phenomenon for which no 
explanation was available. 

Our goal in [3] was thus to develop a deeper understanding of Krause's model 
and of these observed phenomena, by using explicitly the dynamics of the interaction 
topology. To this effect, we introduced a new notion of stability, tailored to such multi- 
agent systems, which provided an explanation for the observed inter-cluster distances 
when the number of agents is large. Furthermore, to understand the asymptotic 
behavior as the number of agents increases, we also studied a model involving a 
continuum of agents. We obtained partial convergence results for this continuum 
model, and proved nontrivial lower bounds on the inter-cluster distances, under some 
conditions. 

Our results in [3] were however incomplete in certain respects. In particular, the 
question of convergence of the continuum model remains open, and some of the results 
involve assumptions that are not easy to check a priori. We see two main reasons for 
these difficulties. First, the system is asymmetric, in the sense that the influence of 
Xj(t) on Xi(t + 1) can be very different from that on Xi(t) on Xj(t + 1), when i and j 



1 The model is sometimes referred to as the Hegselmann-Krausc model. 
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do not have the same number of neighbors. Second, the discrete time nature of the 
system allows, for the continuum model, buildup of an infinite concentration of agents 
with the same opinion, thus breaking the continuity of the agent distribution. 

For the above reasons, we have chosen to analyze here the system (|2.1[) . a con- 
tinuous-time symmetric variant of Krause's model, for which we provide crisper and 
more complete results. One reason is that, thanks to the symmetry, the average value 
n Si x i{t) is preserved, and the average value of a group of agents evolves independent 
of the interactions taking place within the group, unlike Krause's model. In addition, 
when two agent values approach each other, their relative velocity decays to zero, 
preventing the formation of infinite concentration in finite time. The continuous-time 
nature of the system brings up however some new mathematical challenges, related 
for example to the existence and uniqueness of solutions. 

1.1. Outline and contributions. In Section^ we give some basic properties 
of the model (|2.ip . and prove convergence to clusters in which all agents share the 
same value. We then analyze the distance between consecutive clusters building on 
an appropriate notion of stability with respect to perturbing agents, introduced in [3]. 
This analysis leads to a necessary and sufficient condition for stability that is con- 
sistent with the experimentally observed inter-cluster distances, and to a conjecture 
that the probability of convergence to a stable equilibrium tends to one as the num- 
ber of agents increases. In Section [3l we introduce a variant involving a continuum 
of agents, to approximate the model for the case of a finite but large number of 
agents. Under some smoothness assumptions on the initial conditions, we prove the 
existence of a unique solution, convergence to clusters, and nontrivial lower bounds 
on the inter-cluster distances, consistent with the necessary and sufficient for stability 
in the discrete-agent model. Finally, in SectionHJ we explore the relation between the 
two models, and establish that the behavior of the discrete model approaches that of 
the continuum model over finite but arbitrarily long time intervals, provided that the 
number of agents is sufficiently large. 

The results summarized above differ from those those obtained in [3] for Krause's 
model, in three respects: (i) we prove the convergence of the continuum model, in 
contrast to the partial results obtained for Krause's model; (ii) all of our stability and 
approximation results are valid under some simple and easily checkable smoothness 
assumptions on the initial conditions, unlike the corresponding results in [3] which 
require, for example, the distance between the largest and smallest opinions to remain 
larger than 2 at all times; (iii) finally, we settle the problem of existence and uniqueness 
of a solution to our equations, a problem that did not arise for Krause's discrete-time 
model. 

1.2. Related work. Our model (|2.ip is closely related to that treated by Canuto 
et al. [5] who consider multi-dimensional opinions whose evolution is described by 

=J2t(xi(t) - Xj(t)) ( Xi (t) - Xj (t)) , 

3 

or, to a first order discrete time approximation, by Xi(t + 5t) = Xi(t) +Stii(t), where £ 
is a continuouiH nonnegative radially symmetric and decaying function, taking positive 
values only for entries with norm smaller than a certain constant R. Our model is 



2 The continuity assumption appears however unnecessary in the discrete-time case, as was re- 
cently confirmed by one of the authors of [5] in a personal communication. 
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therefore a particular case of their continuous-time model in one dimension with a 
step function for £, except that a step function does does not satisfy their continuity 
assumption. 

The authors of [5] prove convergence of the opinions, in distribution, to clusters 
separated by at least R for both discrete and continuum-time models. Their conver- 
gence proof relies on the decrease of the measured variance of the opinion distribution, 
and is based on an Eulcrian representation that follows the density of agent opinions, 
in contrast to the Lagrangian representation used in this paper, which follows the 
opinion x of each agent. It is interesting to note that despite the difference between 
these two methods for proving convergence, they both appear to fail in the absence of 
symmetry, and cannot be used to prove convergence for the continuum-agent variant 
of Krause's model. 

Finally, the models in this paper are also related to other classes of rendezvous 
methods and opinion dynamics models, as described in [3, 15] and the references 
therein. 

2. Discrete agents. The differential equation (jl.ip usually has no diffcrentiablc 
solutions. Indeed, observe that the right-hand side of the equation can be discontinu- 
ous when the interaction topology changes, which can prevent x from being diffcren- 
tiablc. To avoid this difficulty, we consider functions x : 5R + — > that are solutions 
of the integral version of (jl.ip . namely 




Observe however that for all t at which ii(t) exists, it can be computed using 

2.1. Existence and convergence. Time-switched linear systems are of the 
form x{t) = x(0) + J Q A t x(t) dr, where A t is a piecewise constant function of t. 
They always admit a unique solution provided that the number of switches taking 
place during any finite time interval is finite. Position-switched systems of the form 
x(t) = x(0) + J A x ( t -)x(t) dr may on the other hand admit none or multiple solutions. 
Our model (|2.ip belongs to the latter class, and indeed admits multiple solutions 
for some initial conditions. Observe for example that the two-agent system with 
initial condition x = (—5, \) admits a first solution x(t) = i, and a second solution 
x(t) = xe~ x '. The latter solution satisfies indeed the differential equation (|1.1[) at 
every time except 0, and thus satisfies (|2.1[) . We will see however that such cases are 
exceptional. 

We say that x £ W is a proper initial condition of h2.1\) if: 

(a) There exists a unique x : 5R + — > 5ft™ : t — > x(t) satisfying (|2.1[) . and such that 
x(0) = x. 

(b) The subset of 5R + on which x is not diffcrentiablc is at most countable, and has 
no accumulation points. 

(c) If Xi(t) = Xj(t) holds for some t, then Xi(t') = Xj(t'), for every t' > t. 

We then say that the solution a; is a proper solution of (|2.ip . The proof of the 
following result is sketched in Appendix [5] and a detailed version is available in [4] . 

Theorem 1. Almost all x £ ffl 1 (in the sense of Lebesgue measure) are proper 
initial conditions. 

It follows from condition (c) and from the continuity of proper solutions that if 
Xi(t) > Xj(t) holds for some t, then this inequality holds for all subsequent times. 
For the sake of clarity, we assume thus in the sequel that the components of proper 

4 



initial conditions are sorted, that is, if i > j, then Xi > Xj, which also implies 
that Xi(t) > Xj(t) for all t. Moreover, an explicit computation, which we perform 
in Section [3] for a more complex system, shows that \xi(t) — Xj(t)\ > \x% — Xj\ e~ nt . 
Observe finally that if Xi+i(t*) — Xi(t*) > 1 holds for some t* for a proper solution x, 
then ii+i(t) > and Xi(t) < hold for almost all subsequent t, so that Xi+x(t) —Xi(t) 
remains larg eH than 1. As a consequence, the system can then be decomposed into 
two independent subsystems, consisting of agents 1, . . . , i, and . . . ,n, respectively. 

We now characterize the evolution of the average and variance (sum of squared 
differences from the average) of the opinions. For this purpose, we let F be the set of 
vectors s S 3? n such that for all i, j € {1, . . . ,n}, either Si = Sj, or |sj — Sj\ > 1. We 
refer to vectors in F as equilibria. 

Proposition 1. Let x be a proper solution of The average opinion x(t) = 

— Y^7=i x i(t) * s constant. The sum of squared differences from the average, V(x(t)) = 
(xi(t) — x{t)) 2 , is nonincreasing. Furthermore, with the exception of a countable 
set of times, if Xt £ F (respectively, Xt € F), then the derivative {dV / di)(x{t)) is 
negative (respectively, zero). 

Proof. For all t, except possibly for countably many, 

lx(t) = -^x i (t) = - (x j (t)-x i (t))=0. (2.2) 

dt n ' n f-^ 

Since x(t) is continuous, this implies that x(t) is constant. 

Observe now that, for all t at which x is diffcrcntiable j^V(x{t)) equals 

n n n 

J22(x i (t)-x(t))x i (t)=2^2x i {t)x i (t) = 2^2 Yl Xi(t) { Xj (t) - Xi(t)) , 

i=l i=l i=l j:\xi(t)-Xj{t)\<l 



where we have used the relation (|2.2p twice, and the definition (|2.ip . The right-hand 
side of this equality can be rewritten as 

]T X. l {t){x 3 (t)-X,{t))+ Xj(t)(Xi{t)- Xj (t)), 

liiW-XjW^l \Xj(t)-Xi{t)\<l 

so that 

±V(x(t)) = - J2 {xAt)-Xi{t)f. 

\x t (t)-x } (t)\<l 

The latter expression is negative if x(t) ^ F and zero otherwise. □ 

There are several convergence proofs for the system (|2.ip . We present here a 
simple one, which highlights the importance of the average preservation and symmetry 
properties, and extends nicely to the continuum model. A proof relying on other 
properties and that can be used in the absence of symmetry can be found in [9] . 

Theorem 2. Every proper solution x of 112.1]) converges to a limit x* € F; that 
is, for any i, j, if x* ^ x*, then \x* — x*\ > 1. 

Proof. Observe that by symmetry, the equality 

k 

E E (**•(*) - = 

2=1 j<k, \xi(t)-Xj(t)\<l 



3 The case Xi_|_i(t) — Xi(t) = 1 is more complex. Agents could indeed become "reconnected", as 
in the nonuniqueness example given above because (I2.lt is allowed to fail at countably many times. 
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Figure 2.1. Example of a temporary, "meta- stable," equilibrium. Initially, two clusters are 
formed and do not interact with each other, but they both interact with a small number of agents in 
between. As a result, the distance separating them eventually becomes smaller than 1. The clusters 
then attract each other directly and merge into a single, larger cluster. 



holds for any k and any t. Therefore, it follows from (|2.ip that for all t but possibly 
countably many, 

|5>(*)=E E (Xj(t)- Xi (t)), (2-3) 

i=l i=l j>k, \xi(t)-Xj{i)\<\ 

which is nonnegative because j > k > i implies Xj(t) — Xi(t) > 0. Since yZj—i Xj(t) 
is bounded, this implies that it converges monotonically, for any k. It then follows 
that every Xi(t) converges to a limit x*. We assume that x* k ^ x k+1 an d suppose, 
to obtain a contradiction, that x* k+1 — x* k < 1. Then, since every term Xj(t) — Xi(t) 
on the right-hand side of (|2.3[) is nonnegative, the derivative on the left-hand side 
is asymptotically positive and bounded away from 0, preventing the convergence of 
£<=!*<(*)■ Therefore, x* k+1 - x* k > 1. □ 

2.2. Stable equilibria and inter-cluster distances. By the term clusters, we 
will mean the limiting values to which the agent opinions converge. With some abuse 
of terminology, we also refer to a set of agents whose opinions converge to the same 
value as a cluster. Theorem [5] implies that clusters are separated by at least 1. On 
the other hand, extensive numerical experiments indicate that the distance between 
adjacent clusters is typically significantly larger than one, and if the clusters contain 
the same number of agents, usually close to 2. We believe that this phenomenon can, 
at least partially, be explained by the fact that clusters that are too close to each 
other can be forced to merge by the presence of a small number of agents between 
them, as in Figure 12.11 To formalize this idea we introduce a generalization of the 
system (|2.1[) in which each agent i has a weight Wi , and its opinion evolves according 
to 

Xi(t) = Xi (0) + j Wi(ajj(T)-Xi(T))dT. (2.4) 

J ° j: |x»(t)-x 3 .(t)|<1 
6 



The results of Section 12.11 carry over to the weighted case (the proof is the same) . 
We will refer to the sum of the weights of all agents in a cluster, as its weight. If all 
the agents in a cluster have exactly the same opinion, the cluster behaves as a single 
agent with this particular weiglrjj]. 

Let s € F be an equilibrium vector. Suppose that we add a new agent of weight 5 
and initial opinion xq ■ consider the resulting configuration as an initial condition, and 
let the system evolve according to some solution x(t) (we do not require uniqueness). 
We define A((5, s) as the supremum of \xi{t) — Si\, where the suprcmum is taken over 
all possible initial opinions xq of the perturbing agent, all i, all times t, and all possible 
solutions x(t) of the system (|2.1[) . We say that s is stable if lim^o A(<$, s) = 0. An 
equilibrium is thus unstable if some modification of fixed size can be achieved by 
adding an agent of arbitrarily small weight. This notion of stability is almost the 
same as the one that we introduced for Krausc's model in [2,3]. 

Theorem 3. An equilibrium is stable if and only if for any two clusters A 
and B with weights Wa and Wb, respectively, their distance is greater than d =1 + 

min{VKA , Wb } 
max{ Wa , Wb } * 

Proof. The proof is very similar to the proof of Theorem 2 in [3] . The main idea is 
the following. A perturbing agent can initially be connected to at most two clusters, 
and cannot perturb the equilibrium substantially if it is connected to none or one. If 
it is connected to two clusters A, B, it moves in the direction of their center of mass 
Wa w a A ¥Wb SB ' wnnc t nc two clusters move at a much slower pace, proportional to the 
perturbing agent's weight. We note that, by a simple algebraic calculation, the center 
of mass of two clusters is within unit distance from both clusters if and only if their 
distance is no more than d. 

If the distance between the two clusters is more than d, then the center of mass of 
the two clusters is more than unit distance away from one of the clusters, say from B. 
Therefore, eventually the perturbing agent is no longer connected to B, and rapidly 
joins cluster A, having modified the cluster positions only proportionally to its weight. 
Thus, the equilibrium is stable. 

On the other hand, if the distance between the two clusters is less than d, then 
the center of mass is less than unit distance away from both clusters. We can place 
the perturbing agent at the center of mass. Then, the perturbing agent does not 
move, but keeps attracting the two clusters, until eventually they become connected 
and then rapidly merge. Thus, the equilibrium is not stable. 

If the distance between clusters is exactly equal to d, the center of mass is at 
exactly unit distance from one of the two clusters. Placing a perturbing agent at the 
center of mass results in nonunique solutions. In one of these solutions, the clusters 
start moving towards their center of mass, and the subsequent behavior is the same as 
in the case where the distance between clusters is smaller than d, thus again showing 
instability. Such a solution violates the differential version of (|2.1| only at time t = 
and thus satisfies (|2.ip . □ 

Theorem [3] characterizes stable equilibria in terms of a lower bound on inter- 
cluster distances. It allows for inter-cluster distances at a stable equilibrium that are 
smaller than 2, provided that the clusters have different weights. This is consistent 
with experimental observations for certain initial opinion distributions (see [9] for 
example). On the other hand, for the frequently observed case of clusters with equal 
weights, stability requires inter-cluster distances of at least 2. Thus, this result comes 

4 In the case of non-proper initial conditions leading to multiple solutions, there exists at least 
one solution in which each cluster behaves as a single agent with the corresponding weight. 
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close to a full explanation of the observed inter-cluster distances of about 2.2. Of 
course, there is no guarantee that our system will converge to a stable equilibrium. 
(A trivial example is obtained by initializing the system at an unstable equilibrium.) 
However, we have observed that for a given distribution of initial opinions, and as 
the number of agents increases, we almost always obtain convergence to a stable 
equilibrium. This leads us to the following conjecture. 

CONJECTURE 1 . Suppose that the initial opinions are chosen randomly and inde- 
pendently according to a bounded probability density function with connected support, 
which is also bounded below by a positive number on its support. Then, the probability 
of convergence to a stable equilibrium tends to 1, as the number of agents increases 
to infinity. 

In addition to extensive numerical evidence (see for example [9]), this conjecture 
is supported by the intuitive idea that if the number of agents is sufficiently large, 
convergence to an unstable equilibrium is made impossible by the presence of at least 
one agent connected to the two clusters. It is also supported by results obtained in 
the next sections. A similar conjecture has been made for Krause's model [2,3]. 

3. Agent continuum. To further analyze the properties of (|2.1[) and its behav- 
ior as the number of agents increases, we now consider a variant involving a continuum 
of agents. We use the interval I = [0, 1] to index the agents, and denote by Y the set 
of bounded measurable functions x : I — > 5ft, attributing an opinion x(a) G 5ft to every 
agent in /. As an example, a uniform distribution of opinions is given by x{a) = a. 
We use the function x : I x 5ft + — > 5ft : (a, t) — > x t (a) to describe the collection of all 
opinions at different tiniest We denote by Xt the function in Y obtained by restricting 
x to a certain value of t. For a given initial opinion function xq G Y, we are interested 
in functions x satisfying 

d f 

T Ma)= / {x t {(3)-x t {a))d(3, (3.1) 

where C± C I 2 is defined for any x G Y by 

C i :={{a,f3)el 2 :\x(a)-x((3)\<l}. 

In the sequel, we denote by \x the indicator functions of Cj. 

Note that x , the restriction of x to t = 0, should not be confused with x , 
an arbitrary function in Y intended as an initial condition, but for which they may 
possibly exist none or several corresponding functions x. The existence or uniqueness 
of a solution to (|3.ip is not guaranteed, and there may moreover exist functions 
that satisfy this equation in a weaker sense, without being differentiable in t. For 
this reason, it is more convenient to formally define the model through an integral 
equation. For an initial opinion function xq G Y, we are interested in measurable 
functions x : I x 5ft + — > 5ft : (a, t) — > Xt(a) such that 

x t (a)=x {a)+ [ If (x T (f3)-x T (a))dp)dT (3.2) 

Jo \J P: (a,P)€C XT J 

holds for every t and for every a G /H Similar to the case of discrete agents, one 

5 Note the reversal of notational conventions: the subscript now indicates time rather than an 
agent's index. 

6 A slightly more general definition would require (1 . 1 . 2 1 to be satisfied for almost all a £ /. 
However, this would result in distracting technicalities. 
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can easily prove that for any solution x of (|3.2p . xt := Xt(a) da is constant, and 

Jq (xtipt) — xt) 2 da is nonincreasing in t. 

For the sake of simplicity, we will restrict attention to nondecrcasing (and of- 
ten increasing) opinion functions, and define X as the set of nondecrcasing bounded 
functions x : / — ► 3?. This is no essential loss of generality, because the only quan- 
tities of interest relate to the distribution of opinions; furthermore, monotonicity of 
initial opinion functions can be enforced using a measure-preserving reindexing of the 
agents; finally, monotonicity is preserved by the dynamics under mild conditions. In 
the sequel an element of X will be referred as a nondecreasing function. Further- 
more, if x : / x [0, oo) — * 5ft is such that i ( 6 I for all t, we will also say that x is 
nondecrcasing. 

3.1. Existence and uniqueness of solutions. The existence of a unique so- 
lution to (|3.2p is in general not guaranteed, as there exist initial conditions allow- 
ing for multiple solutions. Consider for example Xo(a) = —1/2 if a G [0,^], and 
xo(a) = 1/2 otherwise. Observe that, similar to our discrete-agent example, xt = xq 
and Xt(a) = xo(a)e~ t are two possible solutions of (|3.2j) . Nevertheless, we will prove 
existence and uniqueness provided that the initial condition, as a function of a, has 
a positive and bounded increase rate; this is equivalent to assuming that the density 
of initial opinions is bounded from above and from below on its support, which is 
connected. 

Our proof of existence and uniqueness is based on the Banach fixed point theorem, 
applied to the operator G that maps measurable functions x : I x [0, ii] — ► 5ft into the 
set of such functions, according to 

(G(x)) t {a) = x (a) + f If {x T (f3) - x T (a)) dp) dr, 

Jt=Q W/3:(a,/3)eC a , T / 



for some fixed t\. Observe indeed that x is a solution of the system (|3.2p if and only 
if xq = x and x = G(x). 

It is convenient to introduce some additional notation. For positive real numbers 
m, M, we call X m C X the set of nondecreasing functions x : I — * 5ft such that 

x((3)-x(a) 
0-a 

holds for every f3 ^ a, and X M C X the set of nondecreasing functions x such that 
for all (3 ^ a, 

x(P) -x(a) < 

We then denote X m l~l X M by X%[, and say that a function x G X is regular if it 
belongs to X^ for some m, M > 0. Let now _Sf be the operator defined on X and 
taking its values in the set of functions from / to 5ft, defined by 

S?(x)(a)= Jx*{a,i)(x(>Y)-x(a))drY. (3.3) 

Observe that ()3.2p can be rewritten as xt(a) = xo(a) + f Q J*C (x T )(a) dr = (G(x))t- 

The proof of existence and uniqueness rests on two important qualitative proper- 
ties of our model. The first, given in Lemma Q] below, establishes that _Sf is Lipschitz 
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continuous on X m . This property will allow us to establish that the operator G is a 
contraction (when t\ is small enough), and to apply Banach's fixed point theorem. 
The second, Lemma [2] below, gives bounds on the rate at which the opinions of dif- 
ferent agents can approach each other. It is instrumental in showing that regularity 
is preserved, allowing us to apply the same argument and extend the solution to 
arbitrarily long time intervals. 

Lemma 1. Let x be a function in X m , where m > 0. The operator J?? is Lipschitz 
continuous at 5 with respect to the \ \ ■ || norm. More precisely, for any yG Y , 

| \J?(x)-J?(y)\L< (2+^)ll*-tflL. 

Proof. Let x G X m , y 6 Y, and 5 — ||5 — yW^. Fix some a £ I, and let 
N x :— {7 : 5(7) — 5(a)| < 1}, N y := {7 : \y(j) — y{a)\ < 1} be the sets of agents 
connected to a, under the configuration defined by 5 and y, respectively. Let also 
N xy = N x n Ny, N x \ y = N x \ N xy , and N y \ x = N y \ N xy . By the definition O 
of we have _Sf (5)(a) = J N ^ {x{j) - 5(a)) d'f and JS?(j/)(a) = j Ny (y(-f) - y{a)) d-y. 
Therefore, 

Jgf(y)(a)-J2f(i)(a) = (y( 7 ) - 5( 7 ) - y(a) + 5(a)) d 7 

+ /iv> -y( a )) d 7-f N , (5( 7 ) ~x{a)) dry. 

It follows from the definition of N x and TVy that |i(7) — 5(a) | < 1 holds for every 
7 G C and |j/( 7 ) — y(ct)\ < 1 holds for every 7 £ N y \ x C A^. This leads to 

|JSf(»)(a) - J2f(£)(a)| < (|y( 7 ) - 5(7)1 + |y(a) - 5(a)|) ^7 + \N x \ y \ + I^WI 
<2\N xy \S+\N x \ y \ + \N y \ x \ 

< 25+\N X \y\ + \Ny\ X \ , 

(3.4) 

where we have used the bound \N xy \ < \I\ = 1 to obtain the last inequality. It remains 
to give bounds on | N x \ y | and | iVj,^ | . 

If 7 £N y \ x , then 7 G N y , and |y( 7 ) — y(cn)| < 1. This implies that 

|s(7) - x( a )\ < 1^(7) - y( 7 )l + |y( 7 ) - ?)( Q, )I + \y( a ) - 5(<*)l < ^ + 1 + & 

Since the same 7 does not belong to N Xi we also have |5( 7 ) — 5(a)| > 1. Thus, for 
every 7 G N y \ x , the opinion 5(7) lies in the set 

[5(a) - 1-25, 5(a) - 1] U [5(a) + l,5(a) + 1+25] , 

which has length at most 45. Since the rate of change of opinions (with respect to the 
index 7) is at least to, we conclude that -AT^a, < 45/m. A similar argument shows 
that jiV^yl < 45/m. The inequality (|3.4p then becomes 

\J*f(y)(a) Sf(x)(a)\ < 25 + 8^ = U + -) \\y - 5^ , 

TO \ TO/ 

which is the desired result. □ 

Lemma 2. Let 5 G Y . Suppose that a, (3 G I , and x(a) < x(/3). Then, 



Jz? (5)(/3) - J5f(5)(a) > - (5(/3) - 5(a)) 
10 



Furthermore, if x £ X m for some m > 0, then 

Sf(x)(J3) - Sf(x)(a) < — (x(B) - x(a)) . 

m 

Proof. Let N a := {7 : \x(j) - x{a)\ < 1} and Np ~ {7 : |x(7) — 5(/3) | < 1} 
be the sets of agents connected to a and f3, respectively. Let now N a p = N a n A^g, 
N a \p = N a \ N a p and JV^\ Q = Np \ N a p. It follows from the definition ([3T3]) of _£f 
that 

JSf (*)(/?) = /^(i( 7 )-ilrf7 + (i( 7 )-5(/3))rf7, 

JSf(5)(a) = J Nag (x(j)-x(a))d 7 + f N ^(x(j)-x(a))d r (3 -°> 

The definitions of the sets Np\ a and N a \p, together with x(/3) > x(a), imply that 
5(7) > x(a) holds for all 7 e Np\ a , and £(7) < i(/3) holds for every 7 e N a \p. Using 
these inequalities and subtracting the two equalities above, we obtain 



Sf{x)(fi) -J2f(5)(o) > / (x(a) - x(fi)) dry + {x{a) - x{p)) dry 

Since |iV aj a| + UVq,^ + |ATg\ Q | = |^Va U iVg| < |J| = 1, we obtain the first part of the 
lemma. 

Let us now assume that x G X m . It follows from (|3.5p and from the inequality 
x([3) > x(a) that 

J8f(s)G3)-JS?(£)(a)< / (i( 7 ) -*(/?)) d 7 - / (i( 7 ) - i(a)) d 7 . (3-6) 



A". 



\0 



which is bounded by UV^wl + UVa\0 • Observe that a^A^^C^a) + 1, £(/3) + 1). 
Since x G X m , we have 

|JV/9\a| < - |«(JV/J\a)| = - (*(/?) - x{a)) . 
m m 

The same bound holds on | N a \p | . The second part of the lemma follows from the 
bound on ([51))) . □ 

The first part of Lemma [2] implies that for any solution x, the difference of the 
opinions of two agents decreases at most exponentially fast. As a consequence, if 
the initial condition of x is an increasing function of a, then xt is also increasing for 
all t. We note that this statement does not necessarily hold for non-regular initial 
conditions initial conditions that are only nondecreasing. 

We can now formally state our existence and uniqueness result, with the rest of 
the proof given in Appendix [B] This result also shows that if the initial condition is 
regular, then the two models given by a differential or integral equation, respectively, 
admit a unique and common solution, which is regular at all times. 

Theorem 4. Suppose that the initial opinion function satisfies xq £ X^f , for 
some m, M > 0. Then the models \3.1\) and 113. 2\) admit a unique and common 
solution x, and x satisfies 

m e-< < Xt{P) - Xt{a) <Me^, (3.7) 
p — a 



for every t and (3 ^= a. 
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3.2. Convergence and fixed points. In this section, we prove that opinions 
converge to clusters separated by at least unit distance, as in the case of discrete 
agents. The proof has some similarities with the one of Theorem [5] It involves three 
partial results, the first of which establishes the convergence of the average value of 
Xt on any interval. Lemma [3] below involves an assumption that xt is nondecreasing. 
By Theorem |4l this is guaranteed if the initial condition is regular. 

Lemma 3. Let x be a nondecreasing solution of the integral equation A3. 2)) . For 
any c £ I, the limit 

lim / xt(a) da 

exists. As a result, the average value (J" fc c Xtia) da)/(c — b) of xt on any positive length 
interval [b, c] converges as t — ► oo. 

Proof. Fix some c £ [0,1] and t\,ti with < t\ < t 2 . The evolution equation 
flESJ) yields 

x t2 (a)da= / x tl {a)da+ / I / / Xx r (ct, (3) (x T (f3) - x T (a)) d(3da I dr, 
'o Jo Jti \Jo Jo J 

(3.8) 

where we have used the Fubini theorem to interchange the integration with respect 
to r and a. We observe that 

Xx T (a,f3) (x T ((3) - x T {a)) df3 da = 0, 

'o Jo 

because of the symmetry property \ Xt (a, (3) = \ Xt (f3, a). Therefore, 

Xx T {»,0) (x T {0) -x T (a)) dfida = / / Xx T (a,P) {x T {P) - x T (a)) d[3da 
la Jo Jo Jc 

(3.9) 

The latter integral is nonnegative, because x T {(3) — x T (a) > whenever a < c < f3. 
Thus, J Q c x t (a)da is a bounded and nondecreasing function of t, hence converges, 
which is the desired result. □ 

Proposition 2. Let x be a solution of the integral equation h3. 2\) such that x± 
is nondecreasing in a for all t. For all a £ I , except possibly for a countable set, the 
limit limt^oo xt{a) exists. 

Proof. Let y(a) = limsup^^ Xt{a). Since, for any t, Xt{a) is a nondecreasing 
function of a, it follows that y{a) is also nondecreasing in a. Let S be the set of all 
a at which y(-) is discontinuous. Since y(-) is nondecreasing, it follows that S is at 
most countable. 

Fix some a S. Suppose, in order to derive a contradiction, that x t (a) does 
not converge to y(a). We can then fix some e > and a sequence of times t n that 
converges to infinity, such that x tn (a) < y(a) — e. In particular, for any 5 > 0, we 
have 

p Ok pa 

/ x tn ([3)d/3< x tn (a) dp = Sx tn (a) < Sy(a) - Se. (3.10) 

J a — 5 J a — 8 

Since y(-) is continuous at a, we can choose S so that y(a — 8) > y{a) — e/3. 
Furthermore, there exists a sequence of times t„ that converges to infinity and such 
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that 



Xr n (a - S) > y{a - 5) - | > y(a) - y . 
At those times, we have 

x Tn (/3)df3> x Tn (a-S)d(3 = 6x Tn {a-6)>8y(a)-—. (3.11) 

a-S Ja-S 3 

However, Eqs. (|3.10[) and (|3.11|l contradict the fact that f _* Xt{a) da converges, thus 
establishing the desired result. □ 

We now characterize the fixed points and the possible limit points of the system. 
Let F C X be the set of nondecreasing functions s such that for every a,j3 G /, either 
s(pt) = s{(3) or \s{a) — s(j3)\ > 1. Similarly, let F be the set of nondecreasing functions 
s such that for almost every pair (a, 0) G I 2 , either s(a) = s(/3) or \s(a) — s((3)\ > 1. 
Finally, we say that s £ I is a /szerf point if the integral equation (|3.2p with initial 
condition s admits a unique solution x t = s for all t. 

Proposition 3. 

(a) Let x be a nondecreasing (in a, for allt) solution of the integral equation LS.ty) . 
and suppose that y(a) = limj_, 00 xt(a), almost everywhere. Then, y £ F. 

(b) If s € F, then s is a fixed point. 

(c) If s is a nondecreasing fixed point, then s G F. 

Proof, (a) We take the limit in Eq. (|3.8[) , as £2 — * 00. Since the left-hand side 
converges, and the integral inside the brackets is nonnegative (by Eq. (|3.9p ). it follows 
that 

liminf / / Xx T (a, f3) (x T (f3) - x T (a)) d[3 da = 0. 
T ^°° Jo Jo 

Using Eq. (|3.9j) and then Fatou's lemma, we obtain 

r-l 



II 



liminf \x T (a, (3) (x T {f3) — x T (a)) df3 da = 0. 



Note that x T (J3)— x T (a) converges to y{(3)—y{a), a.e. If Xy( a > P) = lj then Xx T (&, P) = 
1 for t large enough. This shows that liminfT-^oo \x T (a, P) > Xy( a :P)- We conclude 
that 

1 

Xy (a,f3)(y(f3)-y(a))df3da = 0. 

We integrate this equation over all c G [0, 1], interchange the order of integration, and 
obtain 

1 

Xy(a, (3)(y((3) - y{a)){[3 - a) dp da = 0. 

This implies that for almost all pairs (a, (3), with a < (3, (with respect to the two- 
dimensional Lebesgue measure), we have Xy( a > P) {y(P) ~ y{ a )) = 0, and y(P) > 
y{a) + 1. This is possible only if y G F (the details of this last step are elementary 
and are omitted). 

(b) Suppose that s G F. We have either Xs(&,P) — 0, or s(a) = s(j3). Thus, 
/ Xs{ot, (3)(s((3) — s(a)) d(3 = 0, for all a, and Xt = s for all t is thus a solution of the 
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system. We now prove that this solution is unique. (Recall that uniqueness is part of 
our definition of a fixed point.) 

Since s is bounded and belongs to F, there exists a positive e < 1/2 such that for 
all a, /3 € I, either 5(a) = s((3) or \s(a) — s(/3)\ > 1 + 3e. Let now y be a solution of 
(13. 2p with s as initial condition. Equation (|3.2[) readily implies that \yt{ce) — s(a)\ < e 
for all t £ [0, e] and a £ I. Therefore, for t £ [0, e] there holds \y t {a) - fjt{P)\ < 1 if 
and only if \s(a) — s(/3)\ < 1, and yt is also a solution of the integral problem 

y t (a) = s(a) + f If (y r (/3) - y T {a)) dr) , 

Jt=0 \J p:(a,/3)eC s J 

which unlike (|3.2[) is a linear system since Cg is constant. It can be shown, using 
for example the Lipschitz continuity of the corresponding linear operator, that this 
system admits a unique solution, so that y t — s holds for t £ [0, e]. Repeating this 
reasoning, we obtain y t = s for all t > 0, and s is thus a fixed point. 

(c) Suppose that s is a nondecreasing fixed point. By the definition of a fixed 
point, the function x defined by x t = s for all t is a solution of the integral equation 
(13. 2|) . Since it trivially converges to s and remains nondecreasing, the result follows 
from part (a) of this Proposition. □ 

The following theorem summarizes the convergence results of this subsection. 

Theorem 5. Let x be a solution of the integral equation \3. j|) such that xo is 
regular (or, more generally, such that Xt is nondecreasing for all t). There exists a 
function y £ F such that hmt—xx, Xt (a) = y(a) holds for almost all a. Moreover, the 
set of nondecreasing fixed points contains F and is contained in F. 

3.3. Stability and inter-cluster distances. As in the discrete case, we call 
clusters the discrete opinion values held by a positive measure set of agents at a fixed 
point s. For a cluster A, we denote by Wa, referred to as the weight of the cluster, 
the length of the interval s~ 1 (A). By an abuse of language, we also call a cluster the 
interval s~ 1 (A) of indices of the associated agents. In this section, we show that for 
regular initial conditions, the limit to which the system converges satisfies a condition 
on the inter-cluster distance similar to the one in Theorem [3] From this result, we 
extract a necessary condition for stability of a fixed point. 

Theorem 6. Let xq £ X be an initial opinion function, x the solution of the 
integral equation tS.ty) . and s = limt—.oo x% the fixed point to which x converges. If xa 
is regular, then 

\B-A\>l+ min{WA > WB \ (3.12) 
max.{W A ,W B } 

holds for any two clusters A and B of s. 

Proof. The idea of the proof is to rely on the continuity of Xt (as a function of a) 
at each t to guarantee the presence of perturbing agents between the clusters. Then, 
if (|3.12[) is violated, these perturbing agents will cause a merging of clusters. 

Let A, B be two clusters of s, with A < B. Since s is a fixed point, 1 < B — A. Let 
to = n ^+WB B ^ e tne ^ r cen ter of mass. Condition (|3.12[) is equivalent to requiring 
the center of mass to be at least unit distance away from at least one of the clusters. 
Suppose, to obtain a contradiction, that this condition is not satisfied, that is, that 
to is less than unit distance away from each of the two clusters A and B. 

Since clusters are at least one unit apart, A and B arc necessarily adjacent, and 
since (|3.12j) is violated, B — A < 2. From the monotonicity of xt, there exists some 
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eel such that 

sup{a : lim Xt(a) = A} = c = mf{a : lim xt(a) = B}. 

t — >oo t — >oo 

Moreover, we have the inclusions 

(c-W A ,c) C {a-Aimt^ooXtia) = A} C [c-W A ,c] 
(c,c+W B ) C {a : limt^oo x t (a) = B} C [c,c+W b ]. 



(3.13) 



Let us fix an e > 0. Since Xt{a) converges to s(a) for almost every a, since all 
Xt are nondecreasing, and since clusters are separated by at least one, there exists a 
t' > such that for all t > t' , the following implications are satisfied: 

a < c — Wa — e =>■ x t {a) < A—l, 

a G (c-W A + e,c-e) =4> x t (a) G (A-e,A + e), , 

a G (c+e,c+W B -e) ^ x t {a) G (B-e,B + e), [ ' 

a > c + W B +e x t {a) > B + l. 

We introduce some new notation. To each function x G X, we associate the 
function : — > ( — 1, 1) defined by 



h(q) = / (x(p)-q)d(3, 

Jx- 1 ((q-l,q+l)) 

The value lx(q) represents the derivative of the opinion of an agent whose current 
opinion is q. In particular, the differential equation (|3.1|) can be rewritten as 4zXt(a) = 
l Xt {xt(a)). 

Let us evaluate l Xt (q) for q G [B — 1 + e, A + 1 — e]. (Note that this interval is 
nonempty, because B — A < 2.) Observe first that q — 1 > A — 1 + e > A — 1, because 
£? — A > 1. From the first relation in (|3.14j) and the continuity of x t with respect to 
a, we obtain 

x t (c - W A - e) < A - 1 < q - 1. 

Observe also that q— 1 < A— e. From the second relation in (|3.14|) and the continuity 
of Xt, we obtain 

x t {c- W A + e) > A-e >g - 1. 

A similar argument around q + 1 shows that 

z t (c + W r fe -e)<g+KB + l< x t (c + W B + e). 

Provided that e is sufficiently small, these inequalities and the monotonicity of x t 
imply that 

[c-W A + e,c+W B -e]<Zx- 1 {{q-l,q+l)) C [ c - W A - e,c+W B + e}. 

It also follows from the inclusions (|3 . 14[) that J^_^ A+e (xt(P) — q) d/3 = Wa(A — 
q) + 0(e) 

and j^Y B 6 feG 9 ) ~ ?) ^ = - g) + 0(e). Therefore, 

&*(«) = VF A (A - g) + W B (B -q) + 0(e) = (W A + W B )(m - q) + 0(e), 
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Observe now that since the two clusters do not satisfy condition (|3.12p . their center 
of mass m lies in (B — 1,A + 1). Provided that e is sufficiently small, we have 
me (B — l + e,A+l— e) and therefore, l Xt (B - 1 + e) > and l Xt (A + 1 - e)<0 for 
all t > t'. 

Recall that xo € X^ for some m,M > 0. From Theorem HJ x satisfies the 
differential equation (|3 . 1 1) jjjx t (a) = J£{x t )(a) ~ l Xt (xt(cx)), and also condition (13. 7|) . 
In particular, x t is continuous and increasing with respect to a £ J, for each t. There 
exists therefore a positive length interval J such that Xt'(J) C [B — 1 + e, A + 1 — e]. 
Since ^(.B — 1 + e) > and l Xt (A + 1 — e) < hold for any t > t', and since 
■^xt(a) = l Xt (x t {a)), this implies that x t (J) C [B — 1 + e, A + 1 — e] for all £ > t' . 
Since J has positive length, this contradicts the inclusions (|3.13j) on the convergence 
to the clusters A and B. □ 

We note that the above proof also applies to any solution of (|3.2[) for which xt is 
continuous with respect to a £ I, for all t. 

From Theorem [6j we can deduce a necessary condition for the stability of a fixed 
point, under a classical definition of stability (in contrast to the nonstandard stability 
notion introduced for the discrete-agent system. Let s be a fixed point of (|3.2p . We 
say that s is stable, if for every e > there is a 5 > such that if \\s — Xo\\ 1 < S, 
then \\s — Xt\\ 1 < e for every t and every solution x of the integral equation (|3.2[) 
with xq as initial condition. It can be shown that this classical notion of stability is 
stronger than the stability with respect to the addition of a perturbing agent used 
in Section 12.21 More precisely, if we view the discrete-agent system as a special case 
of the continuum model, stability under the current definition implies stability with 
respect to the definition used in Section |2~2"1 

Corollary 1. Let s be a fixed point of fS.ty) . If s is stable, then for any two 
clusters A and B, 

1+ m^V^M 
maxjW^A, Wb} 

Proof. Suppose that s does not satisfy this condition, and let K be the infimum 
of \\s — s'\\ 1 over all s' satisfying the condition. Clearly, K > 0. For every S > 0, there 
exist M > m > 0, and xq € such that \\s — Sold < Let x be the solution of the 
integral equation (|3.2[) with aio as initial condition, and s' the a. e. -limit of Xt- It follows 
from Theorem |5] that s' satisfies condition (|3.12p . and therefore that \\s — s'Hj > K. 
Using the dominated convergence theorem, we obtain lim^ — >OQ \ \s — — ||^ — ^'lli* 
As a result, lim^oo \ \s — Xt\\ 1 > K > holds for initial conditions xq arbitrarily close 
to s, and s is therefore unstable. □ 

It is possible to prove that the strict inequality version of condition (|3.15|) is also 
necessary for stability. The proof for the case of equality relies on modifying the 
positions of an appropriate set of agents and "creating" some perturbing agents at 
the weighted average of the two clusters. See Chapter 10 of [9] or Theorem 6 in [3] for 
the same proof applied to Krause's model. We conjecture that the strict inequality 
version of condition (|3 . 1 5[) is also sufficient. 

Conjecture 2. A fixed point s of h3.2\) is stable according to the norm \ \ ■ \ \ 1 if 
and only if for any two clusters A, B, 

min{W A ,W B } 
{B ^ l>1+ max{V^,WM' 
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We note that Conjecture [5] is a fairly strong statement. It implies, for example, 
that multiple clusters are indeed possible starting from regular initial conditions, 
which is an open question at present. 

4. Relation between the discrete and continuum- agent models. We now 

formally establish a connection between the discrete-agent and the continuum-agent 
models, and use this connection to argue that the validity of Conjecture [2] implies the 
validity of Conjecture [TJ Toward this purpose, we begin by proving a result on the 
continuity of the opinion evolution with respect to the initial conditions. 

Proposition 4. Let x be the solution of the continuum model i3.2\) for some 
regular initial condition xo G X^ ' . For every e > and T > 0, there exists a 
S > such if y is a solution of the continuum model \3. || ) and \\y$ — ^olloo < $> then 
\\Vt- Xt\L <e, for allte [0,T]. 

Proof. From Theorem|4j xt G X me -t for all t. Lemma[T]then implies that for any 
y t G Y and any t G [0,T], 

\myt)-J?(* t )\L< ( 2 + ^ e *) H^-^Hoo < (2 + ^e T )i|yt-z*IL- (4-i) 

For every a G /, we have 

y t (a) - x t (a) = y Q (a) - x (a) + / (jSf (j/ r )(a) - JSf(ar T ) (a)) dr. 

Jo 

It follows from this relation and from the bound (|4.1[) that 

lltft-*tlL-|ltf-- a! -lloo ^ J ( 2 + ^ eT ) H^ _ ^lloo rfr 
holds for any < s < t < T. This implies that for all t G [0, T], 

lift - ««IL < Htt) - iolL e 4 ( 2+ - eT ) < ||W - SolL e T ( 2+ ^). 

Fix now an e > and take 5 > such that <5e T ( 2+ ™ e ) < e. It follows from the 
inequality above that if | \y — x \ 1^ < S, then \ \y t — XtW^ < e for every £ G [0, T]. □ 

The following result shows that continuum-agent model can be interpreted as 
the limit when n — > oo of the discrete-agent model, on any time interval of finite 
length. To avoid any risk of ambiguity, we use £ to denote discrete vectors in the 
sequel. Moreover, we assume that such vectors are always sorted (i.e., j > i => 
£j > We define the operator G that maps a discrete (nondecreasing) vector to 
a function by G(£)(a) = & if a G l 1 ^-,^), and G(£)(l) = £(n), where n is the 
dimension of the vector £. Let £ be a solution of the discrete-agent model (|2.ip with 
initial condition £(0). One can verify that G(£(£)) is a solution to the continuum- 
agent integral equation for a (|3.2p with G(£(0)) as initial condition. As a result, 
the discrete-agent model can be simulated by the continuum- agent model. The next 
proposition provides a converse, in some sense, over finite-length time intervals. 

Theorem 7. Consider a regular initial opinion function xq, and let (^ n ') n >o be a 
sequence of (nondecreasing) vectors in 5ft™ such that linin^oo | |G(£^ (0)) — xq 1 1 = 0, 
and such that for each n, £' n '(0) is a proper initial condition, admitting a unique 
solution £'"'(£). Then, for every T and every e> 0, there exists n' such that 



G(£<"> (£))-** 
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< € 



holds for all t £ [0,T] and n > n' . 

Proof. The result follows directly from Proposition 2] and from the fact that 
G(£<">(i)) is a solution of fl52J) with the initial condition G(£<™>(0)). □ 

When x is regular, a simple way of building such a sequence (£^(0)) n >o is to 
take £^ (0) = Xq(i/ n) . Theorem[7] implies that the discrete- agent model approximates 
arbitrarily well the continuum model for arbitrarily large periods of time, provided 
that the initial distribution of discrete opinions approximates sufficiently well the 
initial conditions of the continuum model. Now recall that according to Theorem 
[6l and for regular initial conditions, the continuum-agent model converges to a fixed 
point satisfying the inter-cluster distance condition (|3.12|) . The conjunction of these 
two results seems thus to support our Conjecture [TJ that the discrete-agent model 
converges to an equilibrium satisfying this same condition, provided that the number 
of agents is sufficiently large and that their initial opinions approximate some regular 
function. This argument, however, is incomplete because the approximation result in 
Theorem [7] is only valid over finite, not infinite, time intervals. Nevertheless, we will 
now show that this reasoning would be valid, with some exceptions, if Conjecture [2] 
holds. 

Proposition 5. Suppose that xq is regular, and suppose that the limit s of the 
resulting solution x of (3. j[) is stable and its clusters satisfy 

\B-A\>l + min Z A Z B} V ^ 
max{W%VVB} 

Let £(0) €E 5ft™ be a vector whose n entries are randomly and independently selected 
according to a probability density function corresponding to Xq. Then, the clusters 
of the limit of the corresponding solution of (2.1)) satisfy (2.1) ). with probability that 
tends to 1 as n — > oo. 

Proof. Let s = lim^oo x t , which is assumed to be stable and to satisfy (|4.2p . 
Since (|4.2p involves a strict inequality, we see that there exists some K > such that 
the clusters of any fixed point s' that satisfies \\s' — s|| x < K must also satisfy (|4.2j) . 
Furthermore, since s is stable, there exists some e > such that if a solution of the 
integral equation (|3.2|) satisfies \\yt> — s\\ 1 < e for some t', then \\yt — s|| x < K for all 
t > t' . To summarize, if a converging trajectory yt comes within e of s, that trajectory 
can only converge to a fixed point whose clusters satisfy (|4.2p . 

Suppose now that £ (0) is a vector generated at random, as in the statement of the 
proposition, and whose components are reindexed so that they are nondecreasing. It 
follows from Kolmogorov-Smirnov theorems (see [8] for example) and the regularity 
of xq that for any given 8 > 0, the probability of the event ||G(£(0)) — ioll^ < 5 
converges to 1, as n — ► oo. 

Since xt converges to s, a.e., the dominated convergence theorem implies that 
there exists some t' be such that \\xt — s\\ 1 < e/2. Let now £(t) be a solution of (|2.1[) 
for the initial condition £(0), the existence of which is guaranteed with probability 1 
by Theorem[TJ Since G(£(t)) is also solution of the problem (|3.2[) with initial condition 
G(£(0)), Proposition|4]implies that when S is chosen sufficiently small (which happens 
with high probability when n is sufficiently large), we will have ||G(£(i')) — ^t'lloo < - 
e/2, and, consequently, ||G(£(i')) — x# | \ 1 < e/2. Therefore, with probability that 
tends to 1 as n increases, 

\\G(0')) - SW, < ||G(£(0) - st'Hx + \\xt> - % < | + \ = e. 
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It follows that, with probability that tends to 1 as n increases, the limit G(£(i)) is a 
fixed point that satisfies (|4.2p . □ 

We now use Proposition[S]to establish the connection between our two conjectures. 
Suppose that Conjecture holds. Let Xq be a regular initial condition. By Theorem 
[6l the resulting trajectory converges to a fixed point s that satisfies the nonstrict 
inequality (|3. 12jl . We expect that generically the inequality will actually be strict, in 
which case, according to Conjecture[2j s is stable. Therefore, subject to the genericity 
qualification above, Proposition [5] implies the validity of Conjecture [U 

5. Conclusions. We have analyzed a simple continuous-time multi- agent system 
for which the interaction topology depends on the agent states. We worked with the 
explicit dynamics of the interaction topology, which raised a number of difficulties, as 
the resulting system is highly nonlinear and discontinuous. This is in contrast to the 
case of exogenously determined topology dynamics, which result into time-varying 
but linear dynamics. 

After establishing convergence to a set of clusters in which agents share the same 
opinion, we focused on the inter-cluster distances. We proposed an explanation for 
the experimentally observed distances based on a notion of stability that is tailored 
to our context. This also led us to conjecture that the probability of convergence to 
a stable equilibrium (in which certain minimal inter-cluster distances are respected), 
tends to 1 as the number of agents increases. 

We then introduced a variant of the model, involving a continuum of agents. For 
regular initial conditions, we proved the existence and uniqueness of solutions, the 
convergence of the solution to a set of clusters, and a nontrivial bound on the inter- 
cluster distances, of the same form as the necessary and sufficient condition stability 
for the discrete-agent model. Finally, we established a link between the discrete and 
continuum models, and proved that our first conjecture was implied by a seemingly 
simpler conjecture. 

The results presented here are parallel to, but much stronger than those that we 
obtained for Krausc's model of opinion dynamics [3]. Indeed, we have provided here a 
full analysis of the continuum model, under the mild and easily checkable assumption 
of regular initial conditions. 

The tractability of the model in this paper can be attributed to (i) the inherent 
symmetry of the model, and (ii) the fact that it runs in continuous time, although the 
latter aspect also raised nontrivial questions related to the existence and uniqueness 
of solutions. We note however that similar behaviors have also been observed for 
systems without such symmetry. One can therefore wonder whether the symmetry is 
really necessary, or just allows for comparatively simpler proofs. One can similarly 
wonder whether our results admit counterparts in models involving high-dimensional 
opinion vectors, where one can no longer rely on monotonic opinion functions and 
order-preservation results. 

As in our work on Krause's model, our study of the system on a continuum and 
the distances between the resulting clusters uses the fact that the density of agents 
between the clusters that are being formed is positive at any finite time. This however 
implies that, unlike the discrete-agent case, the clusters always remain indirectly 
connected, and it is not clear whether this permanent connection can eventually force 
clusters to merge. In fact, it is an open question whether there exists a regular initial 
condition that leads to multiple clusters, although we strongly suspect this to be 
the case. A simple proof would consist of an example of regular initial conditions 
that admit a closed-form formula for x t . However, this is difficult because of the 
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discontinuous dynamics. The only available examples of this type converge to a single 
cluster, as for example, in the case of any two dimensional distribution of opinions 
with circular symmetry (see [5]). 
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Appendix A. Existence and uniqueness of solutions to the discrete-agent 
equation: Proof of Theorem [T] . 

We sketch here the proof of Theorem [TJ a full version of which is available in [4] . 
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Observe first that if x is the unique solution of the system (|2.ip for a given initial 
condition, then = Xj(i) implies that Xi(t') = Xj(t') holds for all t' > t. Indeed, 
one could otherwise build another solution by switching x\ and Xj after the time t, in 
contradiction with the uniqueness of the solution. Therefore, every initial condition 
Xq satisfying condition (a) automatically satisfies condition (c). 

Let us now fix the number of agents n, and for each graph G on n vertices, with 
edge set E, define X G C 5R ra as the subset in which \x% — Xj\ < 1 if <E E, and 
\xi — Xj\ > 1 if £ E. When restricted to X G , the system (|2.ip becomes the 

linear time invariant differential system 

±i= fe - ^)' ( AJ ) 

which admits a unique solution for any initial condition. This system can be more 
compactly written as x = —L G x 1 where Lq is the Laplacian matrix of the graph G. 

Consider an initial condition i £ S", and suppose that x &X Go for some Go- Let 
xg be the unique solution of x = —L Go x with xg o (0) = x. If this solution always 
remains in Xg , it is necessarily the unique solution of (|2.1[) . Otherwise, let t\ > 
be the first time at which x Go (t) £ dX Go , and set x(t) = XG (t) for all t 6 [0, ti\. By 
the definition of the sets Xq, the point x(t\) also belongs to the boundary of at least 
one other set X Gl , with Gi and Go differing only by one edge We consider here 

the case where £ Eq, $ E\, and Xi(ti) > Xj(ti), but a similar argument 

can be made in the three other possibilities. We also assume that x(ti) belongs to the 
closure of no other set Xq, and that {L Go x{t\))i — {L Go x(t\))j ^ 0. This assumption 
does not always hold, but can be proved to hold for all boundary points that can 
be reached, except for a set that has zero measure (with the respect to the relative 
Lebesgue measure defined on the lower-dinensional boundary). 

Since Xi{t) —Xj(t) < 1 for t just before t\ and Xi(t\) —Xj(t{) = 1, there must hold 
limf[-t 1 (±i(t) — Xj(t)) > 0, and thus —{La x(tx))i + (La x(ti))j > 0, because we have 
assumed that the latter quantity is nonzero. Recall that G\ is obtained from Go by 
removing the edge Since Xi(t±) — Xj(t\) = 1, we have 

-(LG^x)i + {L Gl x)j = -(L Go x)i + (L Go x)j-2(xj(ti)-Xi(ti)) 
= -{LG x)i ; + (L Go x)j + 2 
> 0. 

So, if the solution x can be extended after t\, there must hold Xi{t) — Xj (t) > 1 for all 
t in some positive length open interval starting at ty. This implies that x(t) £ X Gl on 
some (possibly smaller) positive length open interval starting at t±, because x{t\) is at 
a positive distance from all sets Xq other than Xq q and Xq 1 . On this latter interval, 
any solution x must thus satisfy x = —L Gl x. This linear system admits a unique 
solution x Gl for which x Gl (ti) = x{t\). Moreover, the solution remains in X Gl for 
some positive length time period, again because — (La 1 x)i + {L Gl x)j > and because 
x{ti) is at a positive distance from all sets Xq other than X Go and X Gl . If it remains 
in X Gl forever, we extend x by setting x(t) = x Gl (t) on [ti,oo). Otherwise, we 
extend x, as before, on the interval [^1^2], where ti is the first time after t\ at which 
x Gl G dX Gl . In both cases, x is a solution to ()2.ip . on [0, 00) or [0,^] respectively, 
and is unique. Indeed, we have seen that it is the unique solution on [0,<i), that any 
extended solution should then enter X Gl , and that there is a unique solution entering 
X Gl at t\ via x(t\). 
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One can prove that, for almost all initial conditions, this process can be continued 
recursively without encountering any "problematic boundary points," namely, those 
for which (Lcx)i — (Lqx)j = 0, or those incident to more than two sets0 Such a 
recursive construction ends after a finite number of transitions if a solution eventually 
enters and remains forever in a set Xq. In this case, we have proved the existence of 
a unique solution (|2.1|) on 9fJ + , differcntiable everywhere but on a finite set of times. 
Alternatively, the construction may result in an infinite sequence of transition times 
ii, ia, . . .. If this sequence diverges, we have again a unique solution. A problem arises 
only if this sequence converges to some finite time T*, in which case, we could only 
establish existence and uniqueness on [0, T*). The following lemma shows that this 
problematic behavior will not arise, and concludes the proof of Theorem [TJ 

Lemma 4. Suppose that the above recursive construction never encounters prob- 
lematic boundary points (in the sense defined above), and produces an infinite sequence 
of transition times to,t\,.... Then, this sequence diverges, and therefore there exists 
a unique solution x, defined for all t > 0. Moreover, x(t) does not converge to when 
t —¥ oo. 

Proof. Since the sequence t\,t2, ■ ■ ■ of transition times is infinite, a nonempty set 
of agents is involved in an infinite number of transitions, and there exists a time T 
after which every agent involved in a transition will also be involved in a subsequent 
one. Consider now a transition occurring at s\ > T and involving agents i and j. 
We denote by ii(s^) and ii(s^) the limits lim t j Sl ii(t) and lim^ si ii(t) respectively. 
(Note that these limits exist because away from boundary points, the function x is 
continuously differentiable.) 

Suppose without loss of generality that Xi > Xj . If i and j are connected before 
si but not after. The update equation (|2.ip implies that ii(s^) = ii(si) — (Xj(si) — 
Xi(si)). Noting that Xi(si) — Xj(si) = 1, we conclude that that ii(s^) = Xi(s^) + 1. 
Moreover, Xi — Xj must have been increasing just before si, so that ±i(sj~) > ij(si). 
If on the other hand i and j are connected after s\ but not before, then Xj(s^) = 
ij(si) + 1, and since Xi — Xj must have been decreasing just before s±, there holds 
ij(si) > ii(si). In either case, there exits an agent k\ £ {i,j} for which ifc 1 (s^) = 
max{xi(si),Xj(si)} + 1. It follows from si > T that this agent will get involved in 
some other transition at a further time. Call S2 the first such time. 

The definition (|2.1jl of the system implies that in between transitions, < n 

for all agents. Using (|2.1|) again, this implies that ^ 2n 2 for all t at which 

i is not involved in a transition. Therefore, ^^(s^) > irfc 1 (s]'~) — 2n 2 (s2 — Si) = 
x i( s i) + 1 — 2rt 2 (s2 — si). Moreover, by the same argument as above, there exists a hi 
for which ik 2 (sj) = ik x (s^~) + 1 > Xi(s^) + 2 — 2n 2 (s2 — Si). Continuing recursively, 
we can build an infinite sequence of transition times 81,83,... (a subsequence of 
t\, t2, ■ ■ . ), such that for every m, 

ik m i s m) > ii(«r) + m - 2« 2 (s m - Si). 

holds for some agent k m . Since all velocities arc bounded by n, this implies that s m — si 
must diverge as m grows, and therefore that the sequence ti,t2, ■ ■ ■ of transition times 
diverges. □ 

Note that the proof of Lemma [4] provides an explicit bound on the number of 
transitions that can take place during any given time interval. 



7 The authors are pleased to acknowledge discussions with Prof. Eduardo Sontag on this assertion 
and its proof. 
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Appendix B. Existence and Uniqueness solutions to the continuum- 
agent model: Proof of Theorem [4] . 

Let us fix a function £q E X^ , with < m < M . Let us also fix some t\ such 

that 



— < m - IMt < me"* < Me 4t/m < M + —t < 2M, (B.l) 
2 m 

and 

2 + e*<2, (B.2) 

m J 2 

for all t E [0,ti]. We note, for future reference, that t\ can be chosen as a function 
/(m, M), where / is continuous and positive. 

Recall that we defined the operator G that maps measurable functions x : I x 
[0, ii] — > into the set of such functions by 

(G(x)) t (a) = x (a) + f &{x T ){a)dT. 

JQ 

Observe that x is a solution of the integral equation (|3.2p if and only if .To = xq and 
x = G(x). Let P be the set of measurable functions x : Ix [0, t\\ — * 5f : (a, i) — > Tt(a), 
such that .To = To, and such that for all t E [0, ii], we have Tt G X 
detail, 



M+8(Af/m)i 
m-2M( 



m _ 2Mt < x ^l~ x ^ < M + A, (B.3) 
p — a m 



for all t E [0, ti] and all (3 ^ a. In particular, if x £ P, then 

m Xt(0)—Xt(a) „,„ 

— < -^-1 < 2M B.4) 

2 " 0-a ~ V 7 

for all t E [0, ti] and all /3 ^ a. 

We note that P, endowed with the || ■ norm, defined by 

aei,te[o,ti] 

is a complete metric space. We will apply Banach's fixed point theorem to the operator 
G on P. The first step is to show a contraction property of G. 

Lemma 5. The operator G is contracting on P. In particular, \ \G(jj) — G(t)|| oc < 
\\\V- x \\oo f or al1 x,y e P. 

Proof. Let x,y E P. For any t E [0,ii], we have x t E X m / 2 (cf. Eq. (|B.4|) ~). and 
Lemma Q] implies that 

\\^{ yt )-^{ Xt )\\ 00 <(2+-\\\y t -x t \\ 00 . 



m 
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Then, for every a E I, 



\(G(v))t(a) - (G(x)) t {a)\ 



< 



{&{y T )(a)-^{x r )(a)) dr 

f \\^{y T ) - ^{x^dr 
Jo 
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1 ,, 

< 2 Ilf-Slloo' 

where the last inequality follows from Eq. (|B.2[) . □ 

Before applying Banach's fixed point theorem, we also need to verify that G maps 
P into itself. 

Lemma 6. If x e P, then G(x) G P. 

Proof. Suppose that x G P. By definition, G(x)o — xq, and we only need to prove 
that G(x) satisfies condition (|B.3[) . For t G [0, ti] and a < (3, we have 

G(x) t (f3) - G{x) t {a) = x (f3) - x (a) + f (J?(x T )(f3) - JZ(x T )(a)) dr. 

Jo 

It follows from the first part of Lemma [5] and from Eq. (|B.4|) that 

if (x T )(f3) - £?{x T ){u) > - {x T {f3) - x T (a)) > -2M(/3 - a). 
Since xo{(3) — xo(a) > m(/3 — a), for any t G [0,ti], we have 

G(a;) t (/3) - G(x) t (a) > m(/3 - a) - [ 2M(f3 - a) dr = (m - 2Mt){[3 - a), 

Jo 

so that G(x) satisfies the first inequality in (|B.3|1 . 

We now use the second part of Lemma [2] and Eq. (|B.4[) . to obtain 

^(x T )(f3) - J?(x T )(a) < - (x T ((3) - x T (a)) < 8— (fi - a). 

m m 



Since xo((3) — xo(a) < M(j3 — a), for any t G [0, ti], we have 

G(x)t{[3) - G(xUa) < M(f3 - a) + [ 8—(f3 - a)dr =(m+ 8—t) (3 - a). 

n m V m J 



Therefore G(x) also satisfies the second inequality in (|B.3[) . and belongs to P. □ 

By Lemmas [5] and [6j G maps P into itself and is a contraction. It follows, from 

the Banach fixed point theorem, that there exists some a unique x* G P such that 

x* = G(x*). We now show that no other fixed point can be found outside P. 

Lemma 7. If a measurable function x : I x [0,ti) satisfies x = G(x), then it 

satisfies condition \3. 7| ] and, in particular, x G P. 

Proof. Suppose that the function x : I x [0,ti] — ► : (a, t) — > x t (a) satisfies 

x = G(x), that is, xt(a) = xo(a) + J Q ^C(x T )(a) dr for all t and a G I. It follows from 

the first part of Lemma [2] that 

(x t {[3) - x t (a)) - (xo(/3) - x (a)) >— f (x T {(3) - x T (a)) dr 
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holds for a < (3 and < s < t. Together with the fact that x G X m , this implies 
that 

x t ({3) - x t (a) > e~* (x (/3) - x (a)) > me~*(/3 - a), 

which proves the first inequality in (|3.7p and also that Xt G X me -t for all t. Using 
this bound, it follows from the second part of Lemma [2] that 

(x t (P) - x t (a)) - (x ((3) - x (a)) < f* ^e* (x T (/3) - x T (a)) dr 

^ ^l]l( x r{P) ~x T {a))d Tl 

where the last inequality follows from (|B.2[) . Therefore, 

x t (f3) - x t (a) < e 4t / m (i (/3) - so (a)) < Me 4t / m (/3 - a) 

where we have used the assumption that xo G X . This shows the second inequality 
in (|B.3[) and in particular, that x 6 P. U 

We have shown so far that the integral equation (|3 . 2|) has a unique solution x* 
(for t G [0,£i]), which also belongs to P. We argue that it is also the unique solution 
to the differential equation (|3.ip . Since ^£{xX) is bounded for all t, it follows from the 
equality 

x* =x Q + [ 3?(x*)dT, (B.5) 

that is continuous with respect to t under the || ■ || norm, that is, there holds 
lim T ^ t | \x* — x\ 1 1 = for all t. By Lemma [71 x* satisfies condition (|3.7p , and then 
Lemma [1] implies that «5f is a Lipschitz continuous at every x^ . The continuity of x\ 
with respect to t then implies that «if (xl ) also evolves continuously with t. Therefore, 
we can differentiate (|B.5p , to obtain ^Xj(q) = Ji?(xt)(a), for all t and a. The 
function x* is thus a solution to the differential equation (|3.1|) . Finally, since every 
solution of the differential equation is also a solution of the integral equation, which 
admits a unique solution, the solution of the differential equation is also unique. 

To complete the proof of the theorem, it remains to show that the solution x* 
can be extended to all t E 5ft. Let x\ = x^ . It follows from (|3.7|) that x\ £ X^, with 
mi = me" 11 and M\ = Me 4t/, "\ By repeating the argument given for [0, t{\, there 
exists unique x** , defined on / x [t 1; f 2 ], such that 

x** =xi + / ^«*)dr, 

for all £ G [^1,^2], and where ti — t\ = /(mi, Mi). One can easily verify that the 
function obtained by concatenating x* and x** is a (unique) solution of the integral 
and differential equations on [0, £2], and that it satisfies the bound (|3.7|) . Repeating 
this argument, we show the existence of a unique solution on every [0, t n ], with 

t n+1 -t n = f(me- t \Me 4t ^ m ). 

Since this recursion can be written as t n+ i = t„ + g(t„), with g continuous and 
positive on [0, 00), the sequence t n diverges. (To see this, note that if t n < t* for all n, 
then g(t n ) > mino<f<t» g(t) > 0, which proves that t n 00, a contradiction.) This 
completes the proof of Theorem [4] 
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